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We analyze the Rayleigh equation for the collapse of an empty bubble and provide an explanation 
for some recent analytical approximations to the model. We derive the form of the singularity at the 
second boundary point and discuss the convergence of the approximants. We also give a rigorous 
proof of the asymptotic behavior of the coefficients of the power series that are the basis for the 
approximate expressions. 

PACS numbers: 46. 55. dp 

Iii a recent paper Obreschkow et al[l| derived simple and accurate analytical approximations to the solution of the 
Rayleigh equation^ for the collapse of an empty spherical bubble. Their approximants are based on the expansion of 
the solution about the origin of time and can be improved systematically. They showed that those simple analytical 
expressions are suitable for the analysis of experimental cavitation data obtained in microgravity. 

Each approximant is the partial sum of the power series times a function that takes into account the algebraic 
singularity at the other boundary point. The authors also derived an approximate limit of the sequence of partial 
sums in terms of the polylogarithm or Jonquiere's function. To this end, they resorted to a linear fit of the logarithm 
of the expansion coefficients. 

The results obtained by Obreschkow et alQ are partly analytical and partly numerical. In this comment we analyze 
them in a somewhat more rigorous way with the purpose of providing a sound analytical foundation and explanation 
of the main expressions. 

It is sufficient for our purposes to restrict ourselves to the dimensionless Rayleigh equation of motion for a collapsing 
bubbled 

r(t)r(t) + h(t) 2 +e = 0, 

r(0) = 1, r(0) = (1) 



Taking into account the initial conditions and the fact that r(—t) is also a solution we conclude that r(—t) = r(t). 
The equation (JTJ has been written in such a way that the bubble collapses at t = 1; that is to say: r(l) = 0[1|. 
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If we multiply equation ([I]) by r 2 r and integrate with respect to the dimensionlcss time t we obtain 

^f 2 + ^(r 3 -l)=0 (2) 

Note that r(t) = 1 is a solution to Eq. ^ that satisfies the boundary conditions at t = 0. Since this solution does 
not satisfy Eq. ([1} then both equations are not identical. If we solve Eq. ([2]) for dt/dr and integrate between r = 
and r = 1 we obtain the value of 

^3 f 1 r 3 / 2 dr /3^ri 



If, on the other hand, we solve equation J2l) for f 2 and differentiate the result with respect to t we obtain another 

n 

useful equation t L] 

£ 2 

^ = -^4 (4) 
The solution to the Rayleigh equation ^ can be expanded in a Taylor series about the origin as follows: 

oo 

r(t) = j>^' 

g 2 t 2 £ 4 £ 4 19g 6 t 6 59g 8 t 8 4571g 10 t 10 
2 6 180 720 64800 ^ 

Although equations (TTJ), ([2]) and (@| are not identical this series can be obtained from any of them and it converges 
for all < t < 1 as discussed below. 

Taking into account the initial conditions, the behavior of r(t) about t = 1 and the symmetry of the solution 
Obreschkow et al[l| obtained the first and simplest approximant ro(t) = (l — i 2 ) 2 ^ 5 . This expression is considerably 
accurate in a neighborhood of t = because ro(0) = —0.4 is quite close to c\ ~ —0.418. In addition to it, the 
authors found that the error of this expression is smaller than 1% for all t. For this reason they proposed the modified 
power-series approximants 

r„(t) = (l-t 2 ) 2/5 S n (t) 

r, , , 2i 4 - 5£ 2 2 42 - 30£ 2 - 25£ 4 4 

3=0 

where the functions S n (t) are the partial sums for the Taylor expansion of S(t) = r(t)/ro(t) about t = 0. Numerical 
calculation suggests that a,j < and |Oj+i| < \aj\ for all j. Based on these results Obreschkow et al l[ concluded 
that the approximants r n {t) converge monotonically towards r oc (t) as n — > oo. The accuracy of these approximants 
increases with n, but according to Obreschkow et al 1[ r^t) is not identical to r(i) because r(t) in Eq. (H|) and r^t) 
derived from Eq. ([61 do not obey the same asymptotic behavior as t — > 1 . 

Obreschkov et al|l| realized that ln(aj) vs ln(j') is an almost straight line from which they estimated that a,j w 
a±j~ . Based on this approximate relationship they derived the following quite accurate analytical approximation 
to r(t): 

r*(t) = (l-t 2 ) 2/5 [l + ai Lz 2 . 21 (t 2 )] (7) 
where L s (z) = X^^Li z ^ 1 3 s is the polylogarithm or Jonquiere's function. 
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In what follows we will discuss the following points: first, why ro(i) and the approximants of greater order r n (t) 
are so accurate, second, if roo(i) is equivalent to r(t) for all t, and third, why roo(t) is approximately given by Eq. ([7|. 
In order to answer these questions we need the actual behavior of r(t) as t — > 1. 

We can obtain the asymptotic behavior of r(t) as t — > 1 most easily from Eq. (|2J); the result is 

r(t) _ W (1 - Df" _ 5 X 750'/' K(l - t)f" + _ W5 
6 66 

It is worth noting that the leading term 

r a (t) = [im(1 ~ t)]V5 « 1.28371(1 - t) 2 / 5 (9) 
6 

is an exact solution to Eq. ((4]) that does not satisfy the initial conditions. The function (|9]) does not satisfy the other 
two alternative equations ^ and ([2]). 

If we substitute r(t) = ro(t)S(t) into either of the equations ([1]), ([5]) or Q and take the limit t — > l - then we obtain 

= ^--^ — « 0.972867 (10) 
that is consistent with the analytical asymptotic expression ([8]) as follows from 

hm I$L = lim = «^ (11) 

t->i- n,(f) t->i- ro(t) 6 v y 

Obreschkow et alQ already proved that ro(t) is a reasonably good approximation to r(i) in the neighborhood of 
t = 0. Eq. (fTTj) tells us that ro(t) is also quite close to r(t) in the neighborhood of t = 1. For this reason ro(<) is so 
accurate for all t and the approach of Obreschkow et alQ is remarkably successful even when the sequence of partial 
sums S n (t) converges slowly. 

Let us now go into the question whether r n (t) actually gives r(t) when n — > oo. To begin with, note that the 
sequence of partial sums converges for all t < 1 because the singular point closest to the origin is located at t = 1. 
Therefore it is clear that if 5 n (l) converges towards 5(1) as n — > oo then S n (t) converges towards S(t) for all t 
and r n (t) converges towards the actual solution r(t) of the dimcnsionlcss Raylcigh equation. Our numerical analysis 
suggests that S n (l) — > S(l) as n — > oo; compare, for example, 52oo(l) ~ 0.972892 with Eq. (fT0|) . If we accept that 
SooO-) = S(l) then we can easily prove that r^i) satisfies any of the equations H}, © or a s t — > 1. Consider, 
for example, Eq. g]). If we substitute r(t) = r a (t)S(t) then lim t _^i- r(t) 4 r(t) = -24S'(l) 5 /25 = -$ 2 . On the other 
hand lim^- r n (t) 4 r n (t) = -245 , „(l) 5 /25 which proves the point. Therefore, if roo(t) satisfies Eq. ((4| for the most 
unfavorable case t = 1 then it satisfies that equation for all t. 

The approximant Q is quite accurate in the neighborhood of t = 1 because 

lim ^ ( | 2/5 = 2 2 / 5 [1 + ai Li 2 . 2l {l)\ « 1.28363 (12) 

is very close to the exact asymptotic behavior given by Eq. (j9} . In what follows we show how the form of r* (t) emerges 
from the asymptotic behavior of the coefficients <Zj. 

To begin with, note that if f(x) exhibits a branch point at x = xq with exponent a (1 — x/xo) a then the coefficients 
Cj of the Taylor expansion about x = for f(x) behave asymptotically as \cj\ ~ c|aJo| _J j _a_ , where c is a constant. 
Obviously, we are assuming that there is no other singularity closer to the origin or in the vicinity of Xq. For this 
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reason, the coefficients of the original series ([5]) behave approximately as Cj cj~ 7 ^ 5 reflecting the branch-point 
singularity at t — 1 with exponent a = 2/5. Note that the Taylor series about x — for f(x) converges for all 
< x < xq if a > and that the rate of convergence increases with a. The function S(t) exhibits a branch-point 
singularity at t = 1 with exponent a = 6/5 as shown by the asymptotic expansion 

S(t) = ^ = + -^-(1 - t) I 3^(l - 6/5 + 0((1 - 2 ) d3) 

Therefore, the coefficients Oj behave asymptotically as | a-,- 1 ~ aj -6 / 5-1 = aj~ 22 , where a is a constant. This theoret- 
ical result clearly explains the outcome of the linear fitting by which Obreschkow et aljl| obtained the approximant 
r* (t) . The slight discrepancy between the theoretical and numerical exponents is due to the fact that those authors 
fitted all the coefficients a,j and the asymptotic behavior is determined by those of sufficiently large j. If, for example, 
we fit the coefficients aj for 100 < j < 150 then we obtain a much better agreement between theory and numerical 
approximation: |<ij-| « 0.017j~ 2 20 . Obviously, the reason for fitting all the coefficients is the practical purpose of 
obtaining a suitable approximation for all i[l[. In the present case we are mainly interested in explaining the form 
of the approximant ([7]) and for that reason we resort to fitting the coefficients with the largest available orders that 
reflect the asymptotic behavior of r(t) close to t = 1. We also appreciate that the sequence of partial sums S n (t) 
converges for all < t < 1 because 1 + a = 11/5 > 1 and that the rate of convergence of the series ([6]) is greater than 
the one of (0). 

Finally, we want to discuss an alternative power series with much better convergence properties. It is well known 
that in some cases the inverted series exhibits better convergence properties than the original oneQ. The series 



inversion is the basis for the parametric perturbation theory^). In the present case we define the new variable p 

r- 1 2 



ci ¥ 



i 2 z 2 19£ 4 z 3 59£ 6 2 4 4571f 8 z 5 . , 

= z + 1 1 1 \- ... (14) 

3 90 360 32400 v ' 



where z — t , and invert the series to obtain z(p): 

oo 
j=2 



ft 



p 2 e p 3 c 4 p 4 e 7 P 5 e 

= p - — — + ■ L - ^ + — — + +... (15) 

' 3 90 360 10800 v ' 

It follows from the asymptotic expansion 

V6r 5 / 2 VQr 11 / 2 , x 

t = l-- h... (16) 

5£ 22$ V 7 

that z{p) exhibits a singularity of the form (1 — p/po) 5 ^ 2 , where po = 2/£ 2 w 2.3905. Therefore, the coefficients 
bj behave asymptotically as \bj\ ~ ^Ipol - ^ -7 ' 2 , where b is a positive constant. It follows from fitting In ( | foj | ) for 
80 < j < 100 that \bj\ ~ 1.78 x 2.39" J x j" 3 6 , where 3.6 w 7/2 and 2.39 rj 2/^ 2 which confirm the theoretical result. 

Clearly, the coefficients of the inverted series decrease more rapidly than the coefficients of either ([5]) or (|6|). 
Therefore, from a numerical point of view it is convenient to build approximants based on the inverted series. The 
price we have to pay is that the inverted series does not yield r(t) directly, which may not be a serious drawback for 
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some purposes. We can improve the convergence of the inverted series by means of an appropriate summation method 
like the Pade approximants and thus obtain t(p) = \J z(p) which together with r(p) = 1 — ^ yields the parametric 
representation for r(t). 

Summarizing: the most important results of this comment are 

• The particular form and remarkable accuracy of the approximation r»(f) of Obreschkow ct al is now explained 
by the fact that the coefficients cij behave asymptotically as j^ 11 ^ 5 for large j. 

• Present analysis strongly suggests that r 00 (t) is identical to r(t) for all < t < 1 in disagreement with the 
statement of breschkow et al. 

• Alternative approximations to r(t) in terms of the inverted power series exhibit faster convergence than the 
approach of Obreschkow et al, although the inverted power series does not yield r(t) directly 
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